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Abstract. We propose a mathematical derivation of Brinkman's 
force for a cloud of particles immersed in an incompressible viscous 
fluid. Specifically, we consider the Stokes or steady Navier-Stokes 
equations in a bounded domain C M 3 for the velocity field u 
of an incompressible fluid with kinematic viscosity v and density 
1. Brinkman's force consists of a source term favj where j is 
the current density of the particles, and of a friction term Gnvpu 
where p is the number density of particles. These additional terms 
in the motion equation for the fluid arc obtained from the Stokes or 
steady Navier-Stokes equations set in £1 minus the disjoint union 
of N balls of radius e = l/N in the large N limit with no-slip 
boundary condition. The number density p and current density 
j are obtained from the limiting phase space empirical measure 
^EKKiV^i,"*' where Xk is the center of the fc-th ball and Vk 
its instantaneous velocity. This can be seen as a generalization 
of Allaire's result in [Arch. Rational Mech. Analysis 113 (1991) 
209-259] who considered the case of periodically distributed XkS 
with Vk = 0, and our proof is based on slightly simpler though 
similar homogenization arguments. Similar equations are used for 
describing the fluid phase in various models for sprays. 
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1. Introduction 

The subject matter of this paper is the derivation of macroscopic 
models for the dynamics of large systems of solid particles or liquid 
droplets immersed in a viscous fluid (liquid or gas). Specifically, we 
are concerned with the collective effect of the friction force exerted on 
each particle as a result of the viscosity of the fluid together with a no- 
slip condition at the surface of each particle. This type of fluid/solid 
interaction is relevant in several different physical contexts. 

A first example is provided by the sedimentation of solid particles in 
a viscous incompressible fluid (say, a liquid) , typically under the effect 
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of gravity. There is a huge literature on this subject; we shall only 
mention a few, such as [3J, [8], [S]- 

Another example is the case of sprays. Sprays are complex flows 
which consist of a dispersed phase immersed in some viscous fluid. 

Sprays can be described (Cf. [HI [9]) by systems of coupled macro- 
scopic equations (Eulerian-Eulerian modeling) or by the coupling of 
a macroscopic equation and a kinetic equation (Eulerian-Lagrangian 
modeling). 

We do not claim that the results in this paper provide a complete 
derivation of any of these models for sprays, in particular because we do 
not analyze the coupling between the particle and the fluid dynamics. 
Also, our results apply to steady regimes only, for reasons that will be 
discussed below. 

The present work is only aimed at providing a rigorous derivation 
of the Brinkman force created by a cloud of like spherical particles - 
we recall that this force results from the collective effect of the drag 
exerted on the particles by the surrounding fluid. In models for sprays, 
this Brinkman force would typically be responsible for the coupling 
between the motion of the fluid and that of the dispersed phase. 

Our approach of this problem is the homogenization method: we 
more or less follow earlier works such as [7J and pQ which only consid- 
ered periodic distributions of particles. More precisely, the reference [7J 
established the friction term for the Laplace equation in a periodically 
perforated domain with homogeneous Dirichlet boundary condition. 
The case of the Stokes or Navier-Stokes equations was treated in [1] by 
similar arguments. 

The discussion in the present paper differs from pQ in two ways. To 
begin with, only periodic distributions of particles all of which have 
the same velocity (which, by Galilean invariance can be taken as 0) are 
considered in |1] . In the present paper, we consider clouds of particles 
whose phase space empirical measure converges to some smooth phase 
space density. Thus, as long as this (mild) assumption is verified, the 
particles considered here can each have their own instantaneous ve- 
locitjQ. Another difference with [TJ lies in the method of proof, which 
may lead to simplifications here and there. The reference pQ closely fol- 
lowed the argument in [7J by truncating the velocity field in the vicinity 
of each particle, an operation that has the disadvantage of leading to 
velocity fields that fail to satisfy the incompressibility condition. In 



1 We mention also the paper [T3], where an analogous problem is considered 
for the Navier-Stokes equations. This paper (as the references therein) does not 
unfortunately contain any detail about the convergence proof. 



PARTICLES IN A STOKES FLOW 



3 



the present work, the same goal is achieved by removing to the veloc- 
ity field some carefully constructed solenoidal boundary layers so that 
the resulting vector field still satisfies the incompressibility condition. 
Hence the pressure can be integrated out, thereby leading to somewhat 
easier computations and avoiding painful estimates. Yet, our analysis 
borrows a lot from [7] and pQ, especially in the construction of these 
boundary layers. 

We found it convenient to describe the cloud of particles through 
its empirical measure instead of using (marginals of) its iV-particle 
distribution function, as in [3J, [5], [5], [S] — as a matter of fact, most 
of these references assume nearly factorized iV-particle functions, so 
that both viewpoints are essentially equivalent. 

2. Presentation of the model and main results 

2.1. Formal derivation of the model. Consider a system of N iden- 
tical rigid spheres in a viscous incompressible fluid with kinematic vis- 
cosity v and density pf. For simplicity, we assume that the dynamics 
of the spheres is given, and we seek the collective effect on the fluid of 
the drag force on each sphere. We shall make the two following scaling 
assumptions: 

a) the speed of the spheres is assumed to be small enough, so that 
the quasi-static approximation holds for the fluid motion, and 

b) the collective effect of the drag forces exerted on each sphere is 
of the same order of magnitude as the external force field driving the 
fluid. 

First, we outline the quasi-static approximation a). Our starting 
point is the set of Navier-Stokes equations 



where u = u(t, x) G M 3 and p = p(t, x) are respectively the velocity 
and pressure field in the fluid, while r is the radius of the rigid balls 
immersed in the fluid and Xk(t) is the position at time t of the center 
of the k-th ball B Xk m tr . The density of external force per unit of mass 
in the fluid is / = f {i, x) G 1R 3 . 

Notice that, in this model, the effect of solid rotation for each particle 
is neglected — together with the amount of torque particles subject to 
such solid rotations would exert on the fluid. 

Assume that the motion of the spheres occurs at a time scale that 
is long compared to the typical time scale of the external force field /. 
In other words, we postulate the existence of a small parameter r <C 1 



(1) 



d t u + u ■ V x u + V x p 

U(t, -)\dB(x k (t),r) 



v A x u + f , V x -u = 

Xk{t) , 
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such that 

x k (t) = X k (rt) . 

The quasi-static approximation is obtained as follows: defining the slow 
time variable T = rt and u(t,x) = rU(rt,x), the left-hand side of the 
Navier-Stokes equation is rescaled as 

d t u + u- V x u -uA x u = r 2 (d T U + U ■ V X U) - tv A x U . 

Defining 

tF(T,x) = / ( — ,x J , rP(T,x) =p ( — ,x ) , and V k 



3 ~ I 3 ' ~ \— J— J f I 3 -~ I 3 — ' re jrri 

s r / \t J al 

we recast the Navier-Stokes problem above 

r 2 (d T U + U ■ V X U) + rV x .P = t{v A x U + F) , V, • U = , 

f / (^-)U Xfe(T) „=^(T). 

Neglecting all terms of order 0(t 2 ) in (j2J), we arrive at the quasi-static 
Stokes problem 

-v A x U + V X P = F , V* • U = , 
(3) U(T,-)\ dBxk(Thr = V k (T). 

Notice that, in the Stokes problem above, T is only a parameter, so that 
X k (T) and V k (T) can be regarded as independent. In other words, in 
the Stokes problem considered below, it will be legitimate, under the 
quasi-static approximation, to consider X k as a constant and yet to 
allow V k 7^ 0. 

This accounts for item a) above in the derivation of our model; let 
us now discuss item b), namely the collective effect of the drag force 
exerted on the spheres. 

We recall that the drag force exerted on a single sphere of radius r 
immersed in a Stokes fluid with kinematic viscosity u, density p/ is 

Girpfi/rV 

where V is the relative velocity of the sphere — relatively to the speed 
of the fluid at infinity: see [12J §20. 

Hence the collective force field exerted on the fluid by a system of N 
identical such spheres with prescribed dynamics is of the order of 

6irpfi>Nr(V) , 

where (V) is the average relative velocity of the spheres. 

In the sequel, we assume that the parameters v and p/ are of order 
0(1), as well as (V), but we are interested in situations where r C 1 
(small spheres) and N ^> 1 (large number of spheres). In order for 
the collective effect of the immersed spheres to be of the same order 
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as that of the driving external force field, we postulate (without loss of 
generality) that 

(4) Nr ~ Const. > . 

This scaling assumption leads to the mean field approximation listed 
above as b). 

2.2. The quasi-static, mean field limit. Henceforth we use the 
sphere radius as the small parameter governing all limits of interest 
here, and denote it by e > instead of r. Thus we assume that 
N — > oo, e — > and 

(5) Ne = l. 

We further assume that the fluid and the particles considered here are 
enclosed in a domain and denote the volume that is left free for fluid 
motion by 

N 

Q £ = Q\\jB Xk>£ . 

k=l 

In this setting, the Stokes problem for the velocity field u £ and the 
pressure field p e reads 

<6) {- A n^:g; 

Here, the source term g is the ratio of density of external force per unit 
of mass to the kinematic viscosity. This system is supplemented with 
a no-slip boundary conditions for u on the boundary of fl £ : 



(7) 



u \dB XktS = v k , for k = l,..,N, 
u\an = 0, 



where v k is the instantaneous velocity of the (center of mass of the) 
fc-th sphere. 
Denote by 



1 N 



k=i 

the phase space empirical measure of the system of N spheres and by 
(9) p N (x) = / F N (x,v)dv; j N (x) = / F N (x,v)vdv 
its two first moments. 
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It will be convenient to consider the natural extension of u e to Vt 
defined by 



(10) u e (x) 



u e (x) if x G Q £ , 

v k if x G B Xkt£ , k = l,..,N. 



As recalled above, the Stokes' computation of the friction exerted 
on an immersed sphere by the surrounding viscous incompressible fluid 
involves the relative velocity of the sphere to the speed of the fluid at 
infinity. In order to extend Stokes' analysis to the mean field situation 
considered here, we need to assume that the distance between the im- 
mersed particles is large enough compared to their size. Specifically, 
we assume that 

(11) inf \xu — xA > 2r F where r F := e 1 ^ 3 . 

This assumption on the distance between particles is consistent with 
the critical scale for the total number of particles discussed in [10] . The 
assumption (TlTI) allows considering each particle subject to a drag force 
given by Stokes' formula independently of other particles. Obviously, 
we do not know whether (TlTI) is preserved under particle motion, and 
this is why only steady situations are considered here. 

Likewise, we assume for simplicity that the fluid and the particles 
occupy a smooth bounded domain Q C M 3 , and that there is no di- 
rect interaction between the boundary of Q and any of the immersed 
particles: 

(12) inf dist(x fc ,(9fi) > r £ . 

l<k<N 

Theorem 1. Let Q C M 3 be a smooth bounded domain, and consider 
a system of N balls B Xk S for k = 1, . . . , N and e = 1/N included in Q 
and satisfying conditions < f77)) -< f7^) . Assume that the empirical measure 
Fn has uniformly bounded kinetic energy 

sup / / \\v\ 2 F N (x 1 v)dxdv < oo 

N>1 J JnxRZ 

while the macroscopic density and the current converge weakly in the 
sense of measures 

Pn P, 3n as N -> oo 

with p and j continuous on Cl. 

For eachg G(L 2 (fi)) 3 ; letu £ be the unique weak solution in (if 1 (f2 £ )) 3 
of (0), and define u £ as in ( T7^) . Then, u £ converges in (L 2 (Q)) 3 
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to the solution U of 

( -AC/ + Vn = g + 6n(j-pU), 

(13) { V-U = 0, 
[ Olan = 

As a matter of fact, the same techniques as in the proof of Theorem 1 
allow considering the steady Navier-Stokes, instead of Stokes equations. 
The starting point in this case is 

(14) ( ^■Vu £ -uAu £ + Vp £ = g ^ 
^ v ■ u £ — u, 

In writing the system above, we have retained the kinematic viscosity 
v instead of absorbing it in the source term as in the linear, Stokes 
case. Hence, unlike in <^j, g is the density of external force per unit of 
mass (instead of its ratio to the kinematic viscosity). 
The limiting equations in this case are 

( U -VU -v AU + VT1 = g + QTTu(j - pU), 

(15) { V-U = 0, 
[ = 

Let us briefly discuss the uniqueness problem for ffl5|) . By a standard 
energy argument, one finds that, if U\ and U 2 are weak solutions of f[T5"j) . 
they must satisfy 



6vn/ f p\Ur - U 2 \ 2 dx + (U± -U 2 )\\\ 2 



(O) 



< ||VC/i||La(n)||Z7i - U 2 \\ 2 L 4^ 

+ ||^2||L4(Q)||V(L r i - £/2)||L 2 (Q)||t/l - ^211^4(0) 

We first recall (see [TTj p. 9) that 

Il^e||i4(n) < 4||C/ £ || L 2 ( n)||VL' e ||| 2(n) ; 
together with the Poincare inequality, this entails 
(16) ll^ll^(iD<4C7p||Vr7 e ||* a(n) 
where Cp denotes the Poincare constant in the domain Q. Hence 

|2 



v\\V(Ui-U 2 ) llL , m 

I l II Y7fr II - . . Ml XtfTT TT M|2 

2(H) • 



< 2Cp /2 (||Vt/i|| 



Therefore, uniqueness holds for (fl5l) if 

v > 2Cy 2 (||Vf/ 1 |U 2( o ) + ||Vf/ 2 || L2( n)) . 
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But the usual energy estimates for either of the weak solutions U± and 
U 2 shows that 

Hl vc/ jlU 2 (n) < C P (\\g\\ L 2 {Q) + 67n/||j|| £ 2 (n) ) . 
Finally, uniqueness holds for (1151) if 

is 2 >4C 3 p /2 (\\g\\ L 2 m +6nv\\j\\ L 2 (n) ) 
i.e. for v> u = Mhh^n), \\j\\jj>(a),C P ). 

Theorem 2. Under the same assumptions as in TheoremUl and for 
each v > u (\\g\\ L 2^, C P ), consider, for each g G L 2 (fl) and 

each e = 1/N, a solution u £ of the steady Navier-Stokes equations (fl^P 
with the no-slip boundary condition (T7I). Defining its natural extension 
to Q to be u £ as in / TiOj) . one has u £ — > u in L 2 (Q) as e = 1/N — > 0, 
where u is the unique weak solution of (E 



3. Method of proof 

In this section, we present the strategy for the proofs of Theorems 1 
and 2. 

3.1. Introducing correctors. We recall that the weak formulation of 
the Stokes problem flUD-® is 

(17) I Vu £ -VWdx= [ g-Wdx, 

while the weak formulation of the Navier-Stokes problem (THj) . (J7J) is 

(18) v I Vu £ ■ VWdx = u £ ®u £ : VWdx + g- Wdx , 



for each test solenoidal vector field W G (Hq(Q £ )) 3 , i.e. such that 

V- w = o. 

For each w G (V(Q)) 3 such that V • w = 0, we choose test vector 
fields of the form 

W £ = w-B £ [w] 
where B £ [w] G (Hq(Q)) 3 satisfies 

V • B £ [w) = in Q and B e [w) \B Xk , e = w \s x 
Similarly, we approximate the solution by 

U £ = u £ - A e 
where A £ G (Hq(Q)) 3 satisfies 

V • A £ = in Q and A £ \g x e — v k 
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Explicit formulas for the fields A £ and B £ will be given at the end of 
the present section. Notice that, by construction, 

U £ \s Xk ,=W £ k fe , e =0, foraUA: = l,..,iV. 

In addition, the correctors A £ and B £ are chosen so that 

(19) in (H^tl)) 3 , 

(20) A £ ^0 in (H^)) 3 . 
Condition (fT9l) implies that 

W £ ->> w in (^o(fi)) 3 and VW £ ->> Vw in (L 2 (fi)) 9 . 

Moreover, (fT9l) and (1201 imply that 

B 6 H->0 in {L P {Q)) 3 , 

A £ ^0 in (L p (fi)) 3 , 

for each p G [1,6), by the Rellich-Kondrachov compact embedding 
theorem, so that 

W e -»• w in (L p (fi)) 3 for each p e [1, 6) . 

Condition ( 1201) implies that ?7 e and w e behave similarly as e — > 0. In 
the next subsection, we study the asymptotic behavior of U £ , which is 
somewhat simpler to analyze. As we shall see, condition ([2D]) implies 
that 

U e -± U in (H£(Q)) 3 as £ -> 
for both problems (jSJ) and (fT4j) with the boundary condition (j7]). Hence 
C/ £ -»• £/ in (L p (fi)) 3 for 1 < p < 6 

as £ — ► 0. 

3.2. Weak convergence of U £ . Here we show that (some subsequence 
of) U £ converges weakly in (H l {VL)) 3 (assuming (fT9l) and (1201 ). for both 
problems ([6]) and fj!4H with boundary condition (JTj). 

Indeed, for each k — 1, . . . , N, one has U £ \b x lS — 0, so that the weak 
formulation of the Stokes problem becomes : 

l|VC4||i2 (n) = / Vu £ : VU £ dx - I ' VA £ : VU £ dx 
Jn Jn 

g ■ U £ dx - I VA £ : VU £ dx 
u Jn 

< ||fliU 2 (fi)||£4||i,2(n) + ||V^l e ||L2(Q)||VL / £ ||L2(n) 

By the Poincare inequality and (120]) . which entails a uniform bound of 
the form ||V^4 £ || £2(n) < C, we conclude that ||VZ7 e ||i2(n) is bounded. 
Hence there is a subsequence such that U e — ^ U in (H l {£l)) 3 . 
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For the Navier-Stokes problem we have similarly: 

Hl^^lli^n) = v I Vu £ : VU £ dx -u I V A e : VU £ dx 
Jn Jn 

= / g ■ U £ dx - I (u e - Vu £ ) ■ U £ dx - v I VA £ : VU £ dx 
Jn Jn e Jn 

g ■ U £ dx + / U £ ®U £ : VU £ dx -v I VA £ : VU £ dx 



+ / (Ae ® U £ + U £ ® Ae + A £ <8> A) = W £ ofe . 
./fie 

Observe that 

/ U e ®U e :VU e dx= [ U £ -((U £ -V)U £ )dx 
Jn E Jn e 

= \ [ V (U £ \U £ \ 2 ) dx = 
Jn E 

by Green's formula, since U £ \ m = 0. Hence 

Hl^^lli, 2 (n) ^ ||fl'||L2(n)||C4|U2(n) + v||V^4 e ||L2(Q)||V[/ e ||x2(n) 
+ (|| A; ® U £ \\L2(n) + \\U e ® A||l2(q) + HA? ||l 2 (q) 

)||VC/ £ || L 2 (n) 

< Ibllx 2 (n)||t4||i2(n) 

+ (HI VA £ || L 2(^) + 2||A||L 4 (n)||f4||L 4 (n) + II Allien)) || VC/ e ||£2(Q) . 

Applying inequality (|T6l) shows that 
(21) 

(z/-2v / 2Cy 4 ||A||L 4 (n))||Vf/ e || L2{ n ) 

< (Cp||s , IU 2 (fi) + ^||VA e ||L3(n) + ||A|li 4 (n))- 

Recall that || A|U*(fi) ~^ as e — > by (T2T]) . while ||VA e ||i / 2(o) < C by 
(l2"Uj) . Hence the estimate above entails the bound 

\\VU £ \\ L 2 {n) <C. 

3.3. Weak formulations on the whole domain. Next we recast 
the weak formulations (1171) and (1181) in terms of U £ : as we shall see, 
this is somewhat more convenient, at least in taking the mean field 
limit. 

We first discuss the Stokes problem (fTTl) . Observe that 



Vu £ : VWedx = / Vu £ : VW £ dx = g- W £ dx . 
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Expressing u £ in terms of U e , one arrives at 



W £ : VW £ dx + / VA £ : VW £ dx = g- W £ dx 
n Jn Jn 

which, after replacing W £ with its expression in terms of w, leads to 



W £ : Vwdx - / VU £ : VB £ dx + / VA £ : VW e dx = / g ■ W £ dx . 
n Jn Jn Jn 

Transforming the second integral on the right hand side by Green's 
formula, one eventually obtains 

(22) / W £ : Vwdx+ [ U £ -AB £ dx- [ AA £ -W £ dx = [ g-W £ dx. 
Jn Jn Jn Jn 

Under assumptions (fl9l) - (!20l) 

(23) / VU £ : Vwdx -> / W : Vwdx 
Jn Jn 



(24) / g ■ W £ dx -> g ■ wdx 

Jn Jn 

as e — > 0. Thus we are left with computing the limit of 



W e : VBrdx + / VX : VW F dx 



or, equivalent ly, of 



£ . V LJ £ U>Uj * I * »'~ l £ • V £ L 



U E -AB e dx- / A„4 £ -W £ cfe. 



For the Navier-Stokes problem ffl8|) . we follow the same arguments. 
First 

v \ Vm £ ■ VW £ dx — I u £ ®u £ : VW £ dx 

(25) Jue J ^ 

= v Vu £ ■ VW £ dx - / u £ ®u £ : VW £ dx = / g ■ W £ dx 
Jn Jn Jn 

since Vk, k = 1,...,N, are constants and W £ \b x e — 0. Making the 
substitution u £ = U £ + A £ in (1251) . one gets 



+ / g ■ W £ dx + / (A £ ®U £ + U £ ® A £ + A £ ® A £ ) : VW" e cfo 
Jn Jn 
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and, inserting W £ = w — B £ in the equality above, one eventually arrives 
at 

(26) 

v I VU E : Vwdx - v I VU £ : VB £ dx + v I VA : VW £ dx 



U £ ®U £ : VW £ dx + g- W £ dx 
Jn 

+ / (A £ ®U £ + U £ ® A £ + A £ ® A) ■ VW E dx . 
Jn 

Next we pass to the limit as e — > 0; assumptions f|T9|) - fl20|) imply 
that A £ and U £ converge strongly in L 4 (fi), so that 

/ (A £ ®U £ + U £ ® A £ + A £ ® Ae) ■ VW £ dx -> , 



U £ ®U £ : VW £ dx -> U ®U : Vwdx 



5 >0< w e 



Moreover 



! in 
1/ / VC/ e : Vty& -> v VU : Vwdx , 

so that we are left with the task of computing the limit as e — > of 

-1/ / VC/ e : Vi3 £ dx + i/ / VX : VWi. 



At this point, we need to specify how the correctors A £ and B £ are 
constructed. 

3.4. Defining the correctors. Given any smooth function w on B 0s 
and r > s, we designate by \I/ Sir [w] the solution of the following Stokes 
problem : 

" Aty s , r [w} = vn v H , x e B , r \ B , s , 
v-* a>r H = 0, 

**,rM \b , s =W, 

When s = e, r = r F = e 1 / 3 , we define 



(27) 



We denote 7f £ [w] = U ££ i/3[w]), the pressure field associated to ip £ [w] 
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With the function ^eMj we define the corrector B £ as follows: 

N 

(28) B e = Y,^e[w{- + x k )}{x-x k ). 



k=l 



Whenever w is a constant, i.e. w(x) = v, we use the notation 

$ Si r[t)] = \I/ Sir [w] and likewise 4> £ [v] = ip £ [w] . 
With the function <j) £ [v], we define the corrector A £ in the following 



manner: 



N 



(29) 



As = y^ y (j>s[vk\{x - x k ) . 



k=l 



The vector fields A s and B £ so defined are obviously solenoidal ele- 
ments of (Hq(Q)) 3 that verify the conditions 

A e k fc ,«= v k and B £ [w) \s Xk ,= w \b x ^ ■ 
In section [5], we shall prove that A e and B £ verify assumptions (!20|) - (|19|) . 

4. Explicit formulas for the correctors 

The Stokes equations in an annulus can be solved explicitly; in this 
section, we use these explicit formula to express the correctors A £ and 
B £ , and to estimate the quantity 



VU £ : VB £ dx - / VA £ : VW £ dx . 

Occasionally, we will refer to the appendix (section 16. 2p where a few 
standard computations are summarized. 

We start with a formula for — and hence for <j) e [v]. For each 

iGR 3 , denote r = \x\, to — A. Moreover, we denote P^a = (u ■ a)uj 
is the orthogonal projection on the line Rw. 

Whenever 1 < r < R, 



Aa{R)r 2 + 2 (3(R) + 



7 (R) 5(R) 



{I-P u )v 



(30) 
while 



a(R) r 2 + (3{R) + 



l{R) , S(R) 



+ 



&i,r[ v K x ) 



v for x G -Bo,i j 
for x G BZ R . 
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In the formulas above 
3 



(31) a(R) 
while 



9 



s/?:J +0(l/i2*), /3(i?) = — + 0(l/i? 2 ) 



(32) 



j(R) = -f + 0(1/R) , 5(R) = \ + 0(1/R) 



as R — > +oo. 

If one replaces the boundary condition at r = R with the condition 
at infinity 

lim $ = 

\x\^oo 

the solution is 

(33) *U*K*) = i (; + ^) (/-^)« + | 

We denote by n l oo [v] the associated pressure. 

The following relations hold between the pressure fields II i r\v] and 

rii i0 oH : 

(34) 

u ■ V$i,fl[u](x) - ni jfl [u](x)a; = • V$i j00 H(x) - n lj00 (x)a; 



in the limit as R — > +oo. Finally, 



8a{R)r{I-3P u )v+-0[- 



(35) {u-V<S> lj00 [v]{x)-n 1>00 {x)u) = _|(/+3P w )--f (/-3P w )^j 
Using the obvious scaling relation 

e [f](x) = $ ljre / e [v](a;/e) , 
we see that ( 1311 and (l3"Tj) become, for r e [e, r £ ] 



(f>e[v}(x) 



4ai(e) r 2 + 2/? x (£) + 



7i (e) M £ ) 



{I-Pu,)v 



(36) 



with 



ai(e) r 2 + /3 x (e) + 



7i (g) Mf) 

^ y 3 



A(r) {I-P u )v + B(r)P ul v 



(37) ttl (e) = -1 + 0(e 2 / 3 ) , ft(e) = § + 0(£ 4 / 3 ), 



while 



(38) 7i( e ) = -| e + 0(e 5 / 3 ) , ^(e) = i e 3 + 0(e n / 3 ) . 
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Furthermore, for r G [e, r £ ] 

V(f) £ [v] = - (a(r) + b(r)) u <g> {I - P u ) v 



(39) 



with 



+ b(r) 



(I — P u ) V (g UJ + V ■ UJ (I — 3co> £g> Uj) 



(40) a(r) = 6 ( a x r + 



(41) 6(r) = 2a 1 r-£- 3 £. 

We also record the following formulas for scalar products : 
(42) 

V(f) £ [v k ] : V(j) £ [w(x k )} = [(a + b) 2 + b 2 } (v k ■ w(x k ) - P u v k ■ P u w(x k )) 

+ 6b 2 P UJ v k ■ P^w^Xk) , 

and 
(43) 

V<p £ [v k ] : Vw(- + x k ) = - (a + b) uj ■ V(v k ■ w) + buj ■ {v k ■ Vw) 

+ v k ■ u-> (a — 3b)uj ■ (oj ■ Vw) + b v k ■ uj V ■ w . 

In the last formula, we have kept the term V-w, although all the vector 
fields w considered in this work are solenoidal. 



5. Passing to the limit 

First, we prove that the correctors defined in (T2"9"|) converge weakly 
to in Hq(P) in the vanishing e limit. 

5.1. Weak convergence of A £ . Observe that 
(44) 

G(\z\) P z/ \ z \v k ■ P z/ \ z \w(x k )dz = *fv h - w(x k ) / G(r) r 2 dr 

e<\z\<r £ J e<r<r e 

for each function G for which the integral on the right-hand side makes 
sense. Therefore, using ( |36|) and (144)) . we obtain 



(45) Ue[v]\\h {Q) = *f \v\ 2 (jf r 2 (2A 2 + B 2 ) dr + i 
where the last term comes from the integral on Bq^ £ . 
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Since 



r 2 (2A 2 + B 2 ) dr<C (^(ai) 2 (rj - e 1 ) + |ai/9i|(rf - e 5 ) 

+ |a l7 i|(r £ 4 -e 4 ) + (/?i) 2 (r £ 3 -e 3 ) 

+ (|/3i7i| + l«i^i|)(r 2 -e 2 ) 

+ |/3 l7l | log(r e /e) + (71) 2 (r e - e) 



+ I*i7il (e-r )+(*i) 2 (^-i 



1 1 

e r £ 



we obtain from fl37|) -fl38 

(46) 



N 



fc=i 



iV 



2 

L 2 (C) 



L 2 (Q) 



fe=l 



1 * / 

f l^l 2 ( (^ 3 ) + 0(Ne 7 / 3 ) + 0(iV£ 9 / 3 ) + 0{Ne n ^ 1 loge|) 



k=l 







so that A £ — > in L 2 (fi). 

Next, we consider (1421) with u?(xfc) = t>&. Since 

(47) V0 £ H : V0 £ N = [(a + 6) 2 + 6 2 ] \v k \ 2 + [(a + bf + 5 b 2 } (P w v k ) 2 

we obtain 



l|V0 £ 



L 2 (Bo,r £ \So, e ) 



±^K| 2 / [(a + b) 2 + 2b 2 ]r 2 dr 



¥ ( f «? ('I - ^) - 3 7? (r; 1 - e" 1 ) - f 5 2 (rr 5 ' 5 



- 12 ani (r 2 - £ 2 ) + 24 o^i In j - 2 71^ (r £ ~ 3 - e' 

which, together with (!37j) -( l38l) . gives 

(48) \\^^[v k ]\\h [BQtrABQte) <Ce\v k \ 2 . 

Hence, setting r £ = e 1 / 3 , 

N 



\\VA £ \\ 2 L2{n) <CeJ2 l«*| s 



fc=i 



C 




F,v(x,t») M 2 cfc < C". 
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Therefore, there exists a subsequence such that A £ — A in Hq(Q). 
Since *4. e — ► in L 2 (fi), ^4 = and the whole sequence ^4 £ — ^ in 

5.2. Weak convergence of B £ . Next we prove that the sequence of 
correctors given by definition (|28|) converges weakly to in Hq(Q). 

First we estimate ||23 e [w]IU 2 (fi) and 1 1 [iw] 1 1 i 2 (n) ■ In order to do so, 
we consider the solution \I/ S)r [0] of Stokes problem (J27J) with w = 0. 

Using \l/ £i 2£ — ^ e e i/a as test function in problem (|2"T|) with s = e, 
r = e 1//3 , we see that 

(49) ||V* £ie i/3|| L 2 (jBO]reVBoe) < ||V* ej2e ||x2(B 0i2e \B 0ie )- 

According to [T], Lemma 2.2.5, formula 2.2.37, p. 240, the following 
holds for each rj e]0, 1[ and u G H l (B Q i)) 

(50) ||V* fll i[u]||La(B 0l i\Bo,„) < C (\\ Vu\\l*(B 0>1 ) + ^IMU^)^ > 

where the constant C is uniform in 77 and it. Observe that u(x) = 
(j)(2ex) satisfies 

(51) || M || L2(BM) = (2e)- 3 / 2 ||0|| L2( ^ 2£) , 

(52) ||V M || L2(Bo , l) = (2e)- 1 / 2 ||V0|| L2(Bo , 2£) , 
and 

(53) \\V%,2MhHB 0M \B , e ) = (2e) 1 / 2 ||Vvl>,, 1 M|| L2(i?M \B , 1/2 ). 
Using successively Q53J), (H2J, ([50]) and flED, (E2J) , we see that 



|V* £ie i/3[0]|| L2( B Osl/3 \Bo,e) < ( 2 ^) 1/2 l|V%,lM||r»(Bo, 1 \Bb, 1 /») 



<L7((2,) 1 / 2 ||Vn|| L2(Boa) + (l/2) 1 / 2 ( 2£ ) 



1/2 11 n 

IMU 2 (B ,i) 



(54) = C M|V0|U 2(Bo>2e) + (1/2) 1 / 2 (25)- 1 ii^ii^^o,, , 

Assuming that is smooth and 0(0) = implies that 

(55) ||V^[0]|U 2 (b^ 1/3 \Bo, £ ) < Const.e 3 / 2 . 
Since 



N N 



|VB e M|b ( n)<X)ll v & 

[w(a?ib)]|Ug+y^ ||V , 0e[w(-+a?fc)— w(a?fc)]|Ua 



fc=l fe=l 
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it follows from (|48|) with r £ = and v k = w(x k ) and (|55|) that 

||V-B e H|| L 2(Q) < Const. 

Hence there is a subsequence s.t. B £ — ^ B in if^fi). 
On the other hand 

N N 
y^ y 1ps[w{- + X k )]{x - X k ) = ^2(f) £ [w(x k )}(x - X k ) 
k=l k=l 

N 



+ '^2^e[w(- + X k ) - w(x k )](x - X k ) . 



k=l 



By Poincare's inequality and fl55l) . 



N 

^2\\tp £ [w(- + x k ) - w(x k )](x - x k )\\ 2 L 2 

k=l 
N 

< e 2/3 ^2\\Vip £ [w(- + x k ) - w(x k )](x - x k )\\ 2 L2 

k=l 

< CstNe 11/3 . 



Using f Pf6l) with v k = w(x k ) shows that 



N 



B £ = J2^eH-+X k )}^0 



k=l 



in (L 2 (fi)) 3 , so that B = and the whole sequence 

£ £ -0in (H^Q)) 3 . 

Finally, notice that 



N 



|| tp £ [w{- + x k ) - w{x k )]\\ {H i m)3 -> . 

k=l 

5.3. Limit of J n VA £ : W £ . We have 

/ VA £ : VW e cte = / VA £ : VWx - / VA £ : V£ e Mdx 

</n 



(56) 



AT 

^ / V0 £ [w fc ] : V^ e [w(x fc )]rfz 

fc=1 ^e<|2|<r e 
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In view of ( 1201) . we see that 

lim / VA : Vwdx = . 

Moreover, since 

N 

V(J2^[w(- + x k )-w(x k )}) - in (L 2 (Q)) 9 

k=l 

we also have 

N 

/ Vi £ :V(ViK+Xfc)-w(iJd^O. 
^ fc =i 

Next we estimate fl56|) . Recalling ffl3|) . we have 
(57) 



( / 3~(r) r 2 dr J 



r 2 JF( r ) dr = f a\ (r £ 5 - e 5 ) - 6 5? (r e " 5 - £~ 5 ) - f 7l 2 (f - - i) 



V0 £ [f fc ] : V(f) £ [w(x k )]dx = 4irv k - w(x k ) 

£<\z\<r e 

where (according to (jl2|l ) 

(58) ^(r) = ^[(a + 6) 2 + 6 2 ] + 2 6 2 . 
Therefore, 

1 1, 

r £ e' 

-f a l7 i(r 2 -. 2 )-4 7l 5 1 (r £ - 3 -5- 3 ) 

= |e + 0( £ 5 / 3 ). 

Finally, 

(59) / V<p £ [v k ] : V0 e Kar fc )]tfe = 6nev k ■ w{x k ) + 0(e 5/3 ), 

Je<\z\<r £ 

and 

iV 



limY^ / V(p £ [v k ] : V0 £ [w(x fc )]rfx 

£ -°t^A<k|<r s 

= lim ^ ( 6?r £ v k ■ w{x k ) + 0(e 5/3 ) j 
e_> fe=i ^ ' 

= lim (6vr + 0(e 2/3 )) (eiV) v- w(x) F N (x, v) dvdx 



6tt / j(x) • w(x) dx . 
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Therefore, we conclude that 

(60) lim / VA £ : VW £ dx = Qtt j(x) ■ w(x) dx . 

Jn Jn 

5.4. Limit for f n VU e : VB £ . We have 

/ VU £ :VB £ [w}dx = [ VU £ :VCy2^ £ [w{- + x k )-w{x k )])dx 
Jn Jn , , 

(61) 

+ J2 We (' + Xfe ) : VMu>M]dz 

i. i J £<\z\<r F 



k=1 Jc<\z\<r. 

Since 



N 



V(%2i/) e [w(- + x k )-w(x k )]) ^0 
fc=i 

in (L 2 (Q)) 9 , we see that 

N 

/ VU £ : V(J2^[w(- + x k )-w(x k )])dx^0. 
Jn 



k=i 



e<\z\<r £ 



In order to estimate (161]) . we first integrate by parts, denoting by n the 
outward unit normal vector to the sphere dB Xh ^ Te : 

h= W £ (- + x k ) : V<p £ [w(x k )]dz 

J e<\z\<r E 

U £ {- + x k ) ■ A<j) £ [w(x k )]dz 
+ 1 (n- V<f) £ [w{x k )}{- -x k )) ■ U £ dz. 
Next we use the definition of <p £ to compute 

(62) I k = ( n ■ V(j) £ [w(x k )](- - x k ) - 7r e (- - x k ) n J • U £ dz . 

JdB x ., re V / 



At this point, we observe that n = to = jfzf^r (that is, to is "cen- 
tered" on Xk instead on the origin as in section HI). Since V0 e (x) = 
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\ V$i,r e / s (x/e) it follows from that 

h = - [ (w V$i, r ./ e [tw(z fc )](^) - n 1)Jf • C/ £ da; 

£ JdB Xk ,r e \ J 

= - [ fa;-V$i )0O K^)](^)-n 1)OO (2^)a; 
+ 3^) (/ - 3P U ) w(a; fc ) + O ^ )-E/ e dz. 

By (USD, 

h = ~ 2 \ (-!(/ + 3P W ) u;(x fc ) + 3 (/ - 3PJ w(x fe ) 

e 2 \ „ , „ /e 2 



-| (/ - 3PJ w(z fc ) - ) • C/ e dx + O . , 



so that 



4=4 / f-|(/+3Pj'«;(a; fc )+3(/-3Pj W (x fc )yf/ £ rfa;+0 . 

r e JdB Xk ,r e \ J \ r eJ 

Notice that the same result is obtained in (PQ)) by a somewhat different 
procedure. 

At this point, we claim the following strong limits in (if -1 (M 3 )) 3 
that hold for any G G (C;,(R 3 )) 3 — for a proof, see sec. 1 16. ip ) in the 
appendix below: 



N 



(63) '-' 



J>gG(x fc )<W fc ,, e ^p(x)G(x) in (H- x (R 3 )y 
i=i 

J2reG(x k )-uu5 dBxk , re -^f p (x)G(x) in (H-\R 3 )f 



i=l 

Since U £ ^ U converges weakly in (iP(fi)) 3 , we get 

(64) y] / VC/ e : Vip E [w(- + x k )]dz -> - 6tt / w(x)-U p(x)dx . 

k= l Je<\z\<r E JO, 

so that 

(65) lim / VZ7 e : VS e <ix = — 671 / w(x) ■ U p(x)dx. 



5.5. The limit equation. We start from the weak formulations es- 
tablished in section [3] (that is, equations (IT7|) and (118]) ). 
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5.5.1. The Stokes case. In view of the results established in sections 
and 15.41 we pass to the limit in the Stokes problem 

VU ■ Vwdx — g • wdx + Gn I v ■ wfdxdv — w • U pdx 



Therefore, U is the a weak solution of 

( - a u - 6tt (j - p u) + vn = g , 

(66) I V •[/ = (), 

{ U\ dQ = Q. 

Since the problem above has at most one weak solution, the whole 
sequence U £ converges to U in (if 1 (fi)) 3 . 
In addition 

JV 

u £ = U £ + y^ j (f) e [v k ](x - x k ) -»• U 

k=l 

in (L 2 (Q)) 3 , as can be seen from fj46l) . 

This finishes the proof of theorem 1, assuming fl63|) — whose proof 
is deferred to the appendix below 

5.5.2. The Navier-Stokes case. Likewise, for the Navier-Stokes prob- 
lem dHJ), dZD in the limit as e -> 

v I VU ■ Vwdx = / U ®U : Vwdx 



g ■ wdx + 6n v ■ wfdxdv — J w ■ U pdx 

Given p, j and g, there exists u > large enough, so that, for each 
v > v , the problem 

( U-VU-uA U(x) - 6tt (j - pU) + VU = g, 
(67) <^ V-U = 0, 

{ u\ dn = o, 

has a unique weak solution U G (Hq(Q)) 3 . 

Hence the whole sequence U e converges weakly to U in (Hq(Q)) 3 as 
e -> 0. 

As in the Stokes case, (1461) implies that 

N 

U e = U £ + 1 ^2(f) £ [v k ](x -X k ) 
fc=l 

in (L 2 (fi)) 3 . This completes the proof of theorem 2 — assuming again 
that the limits infl63l) hold. 
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6. Appendix 

6.1. Proof of (1631) . We closely follow the method described in [7] in 
the periodic setting. Given G G (Ct(lR 3 )) 3 , we consider two auxiliary 
problems in U&Li B Xk ,r E '■ for k — 1, . . . , N 

-A£ £ = -3G(x k ), 
H \dB Xk , r = r £ G(x k ), 

and 

-A Xe = -G{x k ) -^(& (G{x k ) -00)00 + 2 G(x k )^j+ 3 G(x k ), 
it \dB Xk ,,. E =r £ {G{x k )-uo)uo. 

Next we extend £ e and Xe by in the complement of Ujfc=i Bx h ,r e - Com- 
puting the Laplacian of £ e and Xe in the sense of distributions in the 
whole Euclidean space, we get 



JV JV 

r_ G(xA Finn 



- A£ £ = -3 G(x k ) l Bxk , re + J> G(x k ) 5 9Bi 

k=l k=l 

N 

(68) = -3 N l Bo>re * (Gp N ) + ]T r £ G{x k ) 8 dBxk 



k=l 



N 



-A Xe = -Y,G{x k )l Bxk ^ 

k=l 
N 

(69) + ^r £ (G(x fc ) ■ u) u 6 dBxk , re 

k=l 

N 

= ~ N Wre * ( G Pn) + ^2 r e (G(x k ) ■ U) U S 9Bxkre 

k=l 

+ J! { f ( 6 ( G ( x *) • ^) ^ + 2 G(x fc )) - 3 G(x fc ) | 1 



Bx k ,r e 
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The solutions of the two auxiliary problems above are 



N 



1=1 

Xe(x) = ^ \x-X k \ 2 J (G(x k ) 

while their gradients are given by 

N 

V£ £ (x) = ^2\x-x k \ {u®G{x k ))l Bxk ^, 



^xfc ,re ' 



k=l 
N 



^Xe(x) = E — — (( U ® " UJ ) Ujl B x 

k=l 

N 

+ £ 



TV / 1 | 2 



fc=l 



x ( (G(x k ) -u)I + G(x k ) ® w ) l B;Cfeire . 



with uj 



\x-x k \ ' 

Then, we estimate 



ll&lll^) < Const. r £ 4 J G 2 p N dx = 0(rt), 

||V6||i 2(R 3) < Const. r e 2 J G 2 p N dx = 0(r 2 £ ), 

||Xe||i2(R3) < Const., r 4 J G 2 p N dx = 0(r £ ), 

||Vx e ||£ 2(R3) < Const.tr 2 J G 2 p N dx = 0(r 2 £ ). 

Therefore, C, £ and % £ — > in (H 1 ) 3 , so that both 

(70) A£ £ and A Xe -> in (i/" 1 ) 3 . 

Next, we recall that iV lB , re ~ 1 wea kly in the sense of measures; 
hence 

(71) N i Bors *( GpN )^fp G 
weakly in the sense of measures. Furthermore 

\\Nl Bxk!re *(G PN )\\ L ~<\\G\\ L - 
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so that, by the Rellich compactness theorem, the limit (17TI) holds in 
the strong topology of (H^l) 3 . 

Going back to ( J68l) . we conclude from (1701 that 

N 

Y, r eG{x k )5 dB ^ -> AnpG in (tf^ 1 ) 3 
fc=i 

strongly. 

Next, we apply the same procedure to the second term on the right 
hand side of ( 1691 . 

First, we observe that the last term on that right hand side is 
bounded in L°°(M 3 ) by 11||G||l<», while 
(72) 

J Q {7 ( 6 ( G (^)^V+ 2G (^)) -ZGMyB. h n)dx -> 

for each G (P(IR 3 )) 3 . Applying the Rellich compactness theorem 
again shows that the convergence (172!) holds in the strong topology of 

(ff,; 1 ) 3 

Going back to ( 1691) and using (1701) . ( 17TT) and (1721 shows that 

JV 

5>.(G(a*) • w^a^ - f in (iJ- 1 ) 3 
fc=i 

strongly. 



6.2. Solution of Stokes' problem in an annulus. We first prove 
the explicit formula for in (1301 . (!3T|) and (1521 . by the same method 
as in [12] §20. By symmetry, we seek in the form = 
curl curl(/(r) v ) (where r = Then 

+ -—) (— + -—) f( r ) = Const 
\dr 2 r dr J \dr 2 r dr J 

so that 



f(r) = ar 3 + /5r + 7 + — . 
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Denoting by P^ the orthogonal projection on ui = x/r we arrive at 
formula f[3"Uj) : 



if + j) {I-P.)v-2 f -PuV 



a(R) r 2 + p(R) 



5(R) , >y(R) 



+ 



{I-P u )v 
P„v. 



Because of the boundary conditions, the constants a,/3,7,5 in the for- 
mula above satisfy the following system of equations : 

3 (i? 5 - 1) « + (R 3 -l)/3=i 



5(i? 3 -l)a + 3(R-1)P 
la + lp+± = 5.. 



2' 



-la 



2 ' ' 2 ^ 4 



leading to the estimates (T5TT) and (I52"|) 

3 



a 



8i? 3 



+ R ^ +00 (1/R 4 



7 



1 + O fl _ +00 (l/i2), 



/3 = ^ + 0^ +0O (l/i? 2 
5 = i + O^ +0O (l/i2). 



Next we compute the pressure for the above flow; for simplicity, we 
first write down the following table : 



Vr = uj, 

V{uj-v) = ^, 

UJ ■ P^V = UJ ■ V, 

so that 

V-P^v 



I — UJ(X)U) 



V - P„,V 



Vuj 

r 

uj ■ Vuj = 0, 

V • (a <g> b) = (V • a) b + a ■ V6, 



<E> u; + (uj ■ v) Vui and uj ■ VP W f = 0, 



r 

2 

~u -v, 
r 

v — 3P.,v 



Auj = -uj . 



PARTICLES IN A STOKES FLOW 

To find the pressure, observe that 
AQ ltR [v](x) = 



5a 3(3 3 
20 a - 



r 3 2 r 3 



v- 3 



5 a 3 6 3 



r 3 2 r 3 
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P.U 



so that, up to some unessential additive constant 

3 5a + 3f3 



20 a r 



2r 5 



to ■ P^v 



Letting R — > +oo in the previous expressions leads to formula ([33 
(see also [T2J) : 



(73) $i,oo H («) 

and 

(74) 

Obviously 



Hi, 



2r 2 



W • P, V. 



*i,ooM(i) - [4a(i2) r 2 + 2/3(P)] (J - P w ) v 
2 [a(P) r 2 + /3(P)] P^ + - O ( , 



R 



W'V$yi[l)](l) = w-V$i )0O H(x) 



while 



- 4a(i2)r(2/-P a > + -0 ( - 



n llfl = n ljOC - 20 a(P)r a; -v + \o[ ^ ) . 



P 



as R — > +oo. 

Formula ( 1341) follows as a consequence: 

w • V$i ijR - U 1)R lu = u ■ V$i,oo - rii i00 u; 



8a(P)r(/-3P 1 > + O - 



From (I73|) and (I74p . we arrive at formula (|35|) : 



(75) w • V$i,oo H (x) - ni )0 o u 



(I + 3P. 



P 



(/ - 3P W 



„4' 



Finally, we derive formulas ( 14"2]) . ( 1431) . First rewrite formula ( 139]) for 
V0 £ [f] with r e [e, r £ ] in the form 

(76) V0 e H = -a(r) JV(«) + b(r) [M (v) + v ■ (J - 3 ® w)] 
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where M and iV are two matrix-valued, linear functions of v : 

M(v) = (I — P u ) v®uj — uj®{I — Pj) v, 
N(v) = u®(I-P u )v, 

while 

a(r) = 6(«ir + ^), b{r) = 2 a x r - ^ - 3 -± 
For each v , w G M 3 and each w 6 S 2 

\M(v) : M(w) =JV(v) : N(w) =-M(v) : iV(w) = (v -w-P u v ■ P u w), 
M(v):u}®u] = M(v) : I = N(v) : uj <g> u = N(v) : I = 0, 
lo ® lo : I = u ® u : u ® u = 1, 1:1 = 3. 

Now P2i) and (|43j) follow from (1761) by elementary manipulations in- 
volving the identities recalled above. 
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